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Abstract 



We apply the asymptotic iteration method (AIM) [J. Phys. A: Math. Gen. 36, 11807 (2003)] to 
solve new classes of second-order homogeneous linear differential equation. In particular, solutions are 
found for a general class of eigenvalue problems which includes Schrodinger problems with Coulomb, 
harmonic oscillator, or Poschl- Teller potentials, as well as the special eigenproblems studied recently 
by Bender et al [J. Phys. A: Math. Gen. 34 9835 (2001)] and generalized in the present paper to 
higher dimensions. 
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1. Introduction 

The study of exactly solvable problems in quantum mechanics and the relation between their solu- 
tions represent valuable contributions to mathematical physics. The range of potentials for which 
Schrodinger's equation can be solved exactly has been extended considerably owing to the inves- 
tigations inspired, for example, by super-symmetric quantum mechanics [1], shape invariance [2], 
and the factorization method [3-9]. Coulomb, harmonic oscillator, and Poschl- Teller potentials are 
known examples of exactly solvable problems. It is also known that Coulomb and harmonic oscillator 
potentials are inter-related and that they are members of a family of exactly solvable Schrodinger 
equations [10-12]. In this article, we use the asymptotic iteration method (AIM) [13] to demonstrate 
that Poschl- Teller potentials also belong to this family and that the (Coulomb, harmonic oscillator, 
Poschl- Teller) exact results may be obtained by elementary transformations, as particular cases of 
the following second-order homogeneous linear differential equation: 

,^ (m + 1) 

^ I 1 - 5x^+2 X 



1 - 6x^+2 



(1.1) 



where Q < x < if 6 > and 0<a;<ooif6<0. Here N = -1, 0, 1, ... , and the real 

numbers a , m , and w are to be specified later. 

We should like to make clear at the outset that all the problems discussed in this paper can 
be transformed into special cases of (1.1). The discovery of a common source for the entire set of 
problems and their treatment in arbitrary dimension d by the asymptotic iteration method (AIM) 

are the most significant features of the paper. Our work has two related aspects: the solution of 
differential eqtiations, and the solution of boundary-value problems. In particular we shall show that 
exact solutions of the class of eigenvalue problem of the form 

-y"{x)+x^''+^yix)=Ex''y{x) (iV = -1,0, 1, . . .) (1.2) 

recently studied by Bender et al [14] follow as special cases of the same differential equation (1.1). 
This will allow us to extend the solutions of (1.2) to provide analytic formulae for the solutions of 
the eigenvalue problem 

-y"{x) + fl!}tl±3. + ^2N+2\ y^^^ ^ j^^Ny^^^ (iv = -1, 0, 1, . . .). (1.3) 



Although the problem is posed here in one dimension, it often represents the radial equation for a 
problem in d -dimensions. For example, in the eases of Coulomb and harmonic oscillator problems, 
m(m -|- l)/a;2 might be the angular momentum term in d -dimensions with m = £ + (d — '3>)/2 . 
More generally, we might wish to consider generalizations of these potentials in which m is real 
and non-negative. Indeed, for the Poschl- Teller family of potentials, m is simply a non-negative 
potential parameter. 

The paper is organized as follows. In section 2, we review the asymptotic iteration method. This 
method was recently introduced [13] to construct exact solutions for a wide class of Schrodinger 
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equations and in many cases provided excellent approximate results for nontrivial eigenvalue prob- 
lems [13,15]. In section 3 we investigate the exact solutions of the basic eigenvalue problem (1.1). 
We obtain analytic formulae for the eigenvalues of (1.1) and also explicit expressions for the eigen- 
functions. In section 4, we study the exact solutions of the differential equation (1.3) where the 
connection with the Bender et al class is discussed in detail. In section 6, we developed the exact 
solutions for different classes of Poschl- Teller Potentials. 

2. The Asymptotic Iteration Method (AIM) 

Consider the second-order homogenous linear differential equation 

V" = \o{x)y' + SQ{x)y (2.1) 

for which Ao(.x) ^ and so(.t) arc functions in Coo (a, b) . In order to find a general solution to this 
differential equation we rely on the symmetric structure of the right hand side of (2.1). Differentiating 
(2.1) with respect to a; , we find that 

y'" = Xi{x)y' + s^{x)y (2.2) 

where 

Al = Aq + So + Aq, and si = Sg + sqAo- 
Meanwhile the second derivative of (2.1) yields 

y"" = X2{x)y' + S2{x)y (2.3) 

for which 

A2 = A'l + si + AqAi, and S2 = s'l + sqAi. 
Thus, for (n + 1)*'' and (n + 2)*'' derivative, n = 1, 2, . . . , we have 

= \n-i{x)y' + Sn-i{x)y, and = K{x)y' + Sn{x)y (2.4) 

respectively, where 

A„ = X'n-i + Sn-i + AoA„_i, and s„ = s'^_-^ + soA„-i. (2.5) 
From the ratio of the (n + 2)*'* and (n + 1)*'' derivatives, we have 

We now introduce the 'asymptotic' aspect of the iteration method. For sufficiently large n , if 

^ = 1^ := a, (2.7) 

then (2.6) reduces to 
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which yields the exact solution 



= Ci exp ^ Y^"^*) = ^^^^-i exp /(" + ^o)dt (2.9) 

where Ci is the integration constant, and the right-hand equation follows using (2.6) and (2.8). 
Substituting (2.9) in (2.4) we obtain the first-order differential equation 



y' + ay = Ci exp ^ J {a + Xo )dt 



(2.10) 



which, in turn, yields the general solution to (2.1) as 



y{x) =expi-j adt j C2 + C1 j exp I j {Xoir) + 2a{T))dT j dt 



(2.11) 



3. An exactly solvable class of eigenvalue problem 

In this section we use AIM to investigate the second order homogeneous linear differential equation 

(3.1) 



„ „ , , , (m + 1)\ , , , 
y =2 [ax p{x) y - wp{x)y, 



where a , m and w are real numbers. For special values of a and w , this differential equation can 
immediately be integrated. For example, if w = , then (3.1) yields 

X 2a r* t'p(t')dt' 



/X J I p[z jar 



Furthermore, if p{x) = e' 



-X /w 



2(m+l) rx ^-1 atV»rfi + 



then 



y{x) = e- r "'pW'^* Ci J'' e/' "'f + C2 



We are interested, however, in studying the exact solutions of (3.1) in which p{x) = jti^tjttj N = 
— 1, 0, 1, 2, 3, ... ,. In this case, (3.1) reads 



y" = 2 



rN + l 



[m + 1) \ , wx 



Denote 



Ao(a;) =2 



1 _ X 
ax^+i (m + 1) 



N 



1 _ bx^+2 



, and so{x) = —- 



wx 



N 



l-hxN+^ X J' "'^ l-6a;^+2' 

we may then apply AIM and the asymptotic condition (2.7) yields for n = 0, 1, 2, 3, . . 
• w'^i-l) = n{2a + 2bm +{n + 1)6) for N = -1 



(3.2) 



(3.3) 
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• wZ\0) = 2n(2a + 2bm + {2n + 1)6) for TV = 

• = 3n(2a + 25m + (3n + 1)5) for N = 1 

• w;r(2) = 4n(2a + 25m + (4n + 1)5) for iV = 2 

• w^{3) = 5n{2a + 2bm + (5n + 1)5) for N = 3 

• ... etc. 

Thus, by induction on N , we can easily verify that the ^"(TV) arc given by 



w-{N)=b{N + 2rn{n + p), p=^^I!^±f^, (3.4) 

where N = —1, 0, 1, 2, ... and n = 0, 1, 2, 3, ... . For the exact solutions y„(a;) , we use the generator 
of the exact solutions (2.11), namely 

yn{x) = C2 exp ^ J akdtj , (3.5) 

where n = 0,1,2,... and fc > n is the iteration step mimbcr. Direct computations, using Afe = 
X'i^_i + Sk-i + XoXk-i and Sk = s'j^-i + soXk-i where so and Ao are given by (3.3), imply the 
following: 

• j/o(a;) = 1 since w;^(iV) = where iV = -l,0, 1,2, 

• yi(a;) = -C2(iV + 2)a (1-^x^+2) 

. y2{x) = C2{N + 2fa{a + 1) (l - ?M^x^v+2 + iP+W) ^2iN +2)^^ 

• ... etc. 

Consequently, we arrive at the following general formula for the exact solutions 2/„(x) 

yn{x) = (-l)"C2(Ar + 2)"(a)„ 2-F1 (-n, p + n; a; 6a;^+2) , (3.6) 

where (ct)^ = ^^^^ , a = ^"+^+'^ , and p = . The Gauss hypergeometric function 

2-Fi is defined by 

2F,[-nMc;z) = ±^^^^z-, (3.7) 
fe=0 ^'^''^ 

that is to say, a polynomial of degree n m z . 

4. A Class of exactly solvable eigenvalue problems and their generalization 

An important step in solving eigenvalues problems using AIM [13] is to find a suitable transformation 
that converts the eigenvalue problem under investigation into differential equation of the form (2.1) 
with Aq 7^ . As mentioned earlier (see [13] and [15]), the rate of convergence of AIM is influenced 
by this choice. The asymptotic solution of the eigenvalue problem usually provides a clear indication 
of a suitable transformation to use. For the class of eigenvalue problems of Bender et al 



("^ + ""^y + a'x'^+'^ y{x) = Ex^y{x), < x < 00, 



(4.1) 
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where N = —1, 0, 1, 2, 3, ... , m > —1 , and y{x) satisfying the Dirichlet boundary condition r/(0) = 
, the asymptotic solutions of the differential equation (4.1), for x approaching and oo , suggest 
the following expression of the wave function: 

nT2 



y(a;)=x-+iexp(-^— )/(x), (4.2) 



where f{x) is now to be determined by AIM. Using (4.2) in (4.1), we find that the function f{x) 
has to satisfy the differential equation 



(4.3) 



where w = E — a{2m + N + 3) . This differential equation is a special case of (3.2) with 6 = and 
X G [0, oo) . Thus (3.4) yields 

E = a{2n{N + 2)+2m + N + 3), n = 0,l,2.... (4.4) 

For the general solutions fn{x) ( n = 0, 1, 2, . . . ), it is enough to take the limit in (3.6) as 6 — > . 
Thus the solutions of (4.2) can be written using the limit relation 



in the form 



lim 2Fi{—n, l/b + a; c; zb) = iFi{—n; c; z). (4-5) 

6—^0 

Mx) = {-irC2{N + 2)"(a)„ li^i (-n, a; -^x^'+A , (4.6) 



where a = 2iii±^. 

For = — 1 , (4.1) represents a d -dimensional Schrodinger equation with the Coulomb poten- 
tial, and the eigenvalues are —a? , if m is the angular momentum parameter m = (. + {d — 3)/2 
in d -dimensions. Replacing the parameter E (in this case E becomes a coupling parameter) 
by = Ze^ and = —Enm , we easily recover, from (4.4), the well-known exact eigenvalues 
Enm = — 4[n+m+iy Coulomb problem in d dimensions. Meanwhile, the case N = repre- 

sents the d -dimensional Schrodinger equation for the harmonic oscillator potential and the param- 
eter E itself is now the eigenvalue E = Enm ) where m is again the angular momentum parameter 
m = ^-|- (d — 3)/2. In this case, (4.4) yields the well-known exact eigenvalues Enm — a{4:n + 2ni + 3) 
for this harmonic oscillator problem in d dimensions. The parameter m may be more generally 
any non- negative real number. For N > we emphasize that the eigenvalues of the boundary- value 
problem (4.1) do not correspond to those of a conventional Schrodinger equation. However, the zero 
eigenvalues may be the ground-state energies corresponding to an effective potential 



F(a.) = ^(^ + ^) +aV^+^-ffx^, ig>0). 



Indeed our method gives us the conditions on m , a , and g under which such zero eigenvalues 
might occur. For example, for N = 2 , we have for 

-^ + ^^2 — - + ax^-gx^)y{x)=Ey{x), < x < oo (4.7) 
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the energy is zero E = if 



g = a{8n + 2m + 5), n = 0, 1, 2, . . . . (4.8) 



5. Poschl- Teller Potentials 

In this section wo show how the asymptotic iteration method can be used to generate the exact 
solutions of the different famihes of Poschl- Teller potentials [16-18] 

(0<fcu< |, a,p>0) (5.1) 

(0 < fcu < oo,a > /3). (5.2) 

Our approach shows that the exact solutions of Poschl- Teller potentials follow directly from solu- 
tions of the differential equation (3.2) through elementary transformations. The main idea is to use 
trigonometric or hyperbolic mappings that make the potentials (5.1) and (5.2) rational and then 
make the corresponding Schrodinger equation suitable to be analyzed with AIM. 



and 



Vuiu) = fc2 



aja + l) ^ /?(/? + 1) 
cos2(A:u) sin^(fcM) 

Pif3-1) aia + 1) 



sinh^(fcw) cosh^(fcw) 



5.1 The first class of Poschl- Teller potentials 

We consider the one-dimensional Schrodinger equation of a quantum particle trapped by the 
Poschl- Teller potential Vi{u) 

-^ + V,{u)y = £y (5.3) 

where Vi{u) is given by (5.1). The potential Vi{u) is closely related to several other potentials 
which are widely used in molecular and solid state physics like the symmetric Poschl- Teller potential 
a = /3 > , and the Scarf potential — ^ < a < . The substitution x = ku , < x < ^ , yields 

The potential Vi{x) = a{a + l)/cos^(a;) -|- /?(/? -|- 1)/ sin^(a;) is periodic. However, each period is 
separated from the next by an infinite potential barrier so it can be studied within one period, say 
< X < ^ . Clearly, (5.4) is a smooth approximation for a, /? ^ 0^ of the infinite square-well 
potentials over the interval [0, -j] . In order to apply the asymptotic iteration method, discussed 
earlier, to solve the eigenvalue problem (5.4), we have to transform it into a differential equation of 
the form (2.1) with suitable Aq 7^ . The asymptotic solutions of (5.4) as x approaches 0+ and 
^~ suggest for y{x) the following expression 

y{x) = cos"+i(x) sm'^+^{x)f{x) (5.5) 

where f{x) is found by means of AIM. After substituting (5.5) into (5.4), we obtain the following 
differential equation for / 



/" = 2{{a + 1) tan(a;) -{l3+l) cot{x)) f - wf 



(5.6) 
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where w = E — {a + P + 2)"^ . Further, if we use the substitution t = cos(x), < t < 1 , we arrive 
at the following equation 

which we may compare with (3.2) for 6=1, a = /3+|, m = a, and N = . Thus from (3.4) we 
have w = 4n(a + (3 + n + 2) where n = 0, 1, 2... . Finally, the eigenvalues for Eq.(5.4) are given by 

En = {a + f3 + 2 + 2nf, n = 0,1,2, ... . (5.8) 

Directly solutions /„(i) for the eigenvalue problem (5.8) can now be obtained from (3.6) with the 
substitution 6=1 and N = . Thus, we have 

fn{t) = (-2)"(:72 {(T)n iFx{-n, p + n,a;t^), t = cos{x) (5.9) 

where p — a + P + 2 and cr = a + | . Consequently, the wave functions for the first Poschl- Teller 
Potential (5.4) are 

yn{x) = (-2)"C2(a + |)„ cos"+^(a;) smf^+\x) 2Fi{-n, a + f3 + 2 + n,a +1; cos\x)) 

2 ^ ^ 2 (5.10) 

= (-2)"(72 n! cos«+i(a;)sin''+^(x) Pi"^^'^^^\l + 2cos^ x) 

using the relation [19] between the hypcrgeometric function and Jacobi polynomial Pn°''^\x) of 
degree n . The constant C2 in (5.10) still to be determined by normalization. 



5.2. The second class of Poschl- Teller potentials 



We now consider the Schrodinger equation 



-^ + Vii{u)y = £y, < w < 00 (5.11) 



where Vii(u) is given by (5.2). With the substitution, x = ku , we have 

d^y (I3{(3-1) a{a+l) 



dx"^ Vsinh^(a;) cosh^(a;) 



y = Ey, < x < 00 (5.12) 



We may assume a > /3 , for if a < p , we can change a — > — a — 1 , because the equation remains 
unchanged under a — > —a — 1 and /3 — /3 + 1 . Because of the asymptotic solutions for y as x 
approaches 0+ and 00 , we may assume the exact form of the wavefunction y{x) to be 

y{x) = cosh-'*(x) sinh''(a;)/(a;), (5.13) 

where the function f{x) is to be fotmd with AIM. After substituting (5.13) into (5.12), we obtain 
the following differential equation for f{x) 

f" = 2 (atanh(a;) - /?coth(a;)) /' - wf, (5.14) 
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where w = E + {a — . Further, by means of the substitution t = sinh(a:) , we can easily obtain 

where the prime refers to the derivative with respect to the variable t . This differential equation 
can be compared with (3.2) with b — —1 , a = a— ^ , m = (i—l and = . Thus, using (3.4), we 
find the following expression for w : w = 4n(a — /? — n) , where n = 0, 1, 2... . Finally, the eigenvalues 
for (5.12) are given expHcitly by the expression 

E„ = -{a- P-2nf, n = 0, 1, 2, . . . < (a - /3)/2. (5.16) 

Consequently, the expected maximum number of quanta is 

n^a. = [(a - /3)/2] (5.17) 

where [(a — /?)/2] stands for the closest integer to {a — (3)/2 that is smaller than (a — /3)/2. 

With the substitution b = —1 and = , wo find using (3.6) that the exact solutions of fn{y) 
which satisfy (5.14) for n = 0, 1, 2, . . . < (a — /?)/2 are 

/„(y) = (-2)"C2 (/3+J)„2Fi(-n,/3-a + n;/3+ J;-t2) ^ ^ 

^ ^ , (5.18) 

= (-2)"C2 n! P}r~^'^~^\l + 2t'^) 

where in the last line we have expressed the hypergeometric function in terms of the Jacobi poly- 
nomial in order to make easier the comparison with the results found in the literature [20]. Up to a 
normalization constant the wave function y{x) of (5.2) reads 

2/„(:c) = (-2)"C2(/3+i)„ cosh-"(x)sinh''(x) 2Fi{-n, P - a + n; (3 + - smh\x)). (5.19) 



6. Conclusion 

In this paper we have applied the asymptotic iteration method to obtain the exact solutions for 
an interesting class of differential equations (3.2). We have shown that the exact solutions for the 
class of problems studied by Bender et al are obtained by an application of AIM to a transformation 
of (3.2). This allows us to obtain the exact solutions for an extended class of eigenvalue problems 
(4.1) with m > —1 . Complete solutions for the well known coulomb and harmonic oscillator po- 
tentials follow directly by setting N = —1,0 respectively in (4.1). By means of trigonometric and 
hyperbolic coordinate transformations of t , in terms of which the first and second kind of Poschl- 
Teller potentials become rational functions of t , AIM also provides a direct way of generating the 
exact solutions to the Schrodinger eigenvalue problem generated by these potentials. 
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